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In a recent Letter ||l[, Cuccoli et al presented a new 
theoretical approach to the understanding of the two 
dimensional (2D) quantum Heisenberg antiferromagnet 
(QHA). The approach is characterized by the feasibility 
of the separation of the quantum fluctuation from the 
thermal one. One of their main predictions is that f for 
the QHA is given by that of the classical counterpart, 



I.e. 



m^iciitci), t,i = t/0\t), 



(1) 



where t = kT/J for S — 1/2 and 9{t) can be determined 
from their theory using their Eqs.(6)-(8), at least numer- 
ically. 

The authors concluded that their theory "appears to 
explain all the experimental data for different values of 
S without any fitting parameters" . Consequently, they 
questioned the validity of the key assumption of the con- 
ventional theoretical approach [||. This conclusion, how- 
ever, is surprising because the Gaussian approximation 
they used to handle the quantum fluctuation is supposed 
to be accurate only when the effect of the fluctuation is 
weak, i.e., for sufficiently large value of S and high tem- 
perature. Since the experimental data have rather large 
statistical fluctuation and previously available MC data 
limited to ^ < 30 are not completely free of systematic er- 
rors, much more accurate data including larger ^ appear 
to be crucial to check the validity of their theory. 

In this Comment we present high precision MC data 
of the second moment correlation length for S — 1/2 
QHA on L X L lattices, up to ^ = 95.7(3). Our data are 
obtained using a powerful new quantum MC method g] 
which completely eliminates the systematic error coming 
from finite Suzuki- Trotter number. We carefully moni- 
tored the finite size effect in our data by repeating mea- 
surements on varying lattice size from L— 20 to 1000 
and found that it becomes smaller than the typical sta- 
tistical error of 0.3 percent or better under the condition 
i/C > 7. 

The classical correlation length £,ci{tci) can be deter- 
mined from extensive MC data already available for the 
2D classical Heisenberg ferromagnet [|| much more pre- 
cisely than using the data given by Cuccoh et al. [Q, 
because an extrapolation oi ^ci{tci) for smaller td is nec- 
essary. Considering ^c/ over the range 34 < S^ci ^ 788 
(corresponding to 0.20 ^ t < 0.29), we find a good 
fit assuming the functional form from the two-loop or- 
der of the perturbation theory for the classical model, 
i.e., Cc/(td) = Q ic/ e°/*=' [1 + d t^ + . . .] , with Q 



2.300 X 10^ 



6.378, ci ~ -0.855, and x /Ndf (the 



X^ value per degree of freedom) ~ 0.4. For the ^c/ over 
2 < Ccz < 34 (corresponding to 0.29 < t < 0.65), on the 
other hand, it turns out that an additional correction of 
the type 02^^; is necessary for a good fit. The estimated 
fitting parameters in this case are: Cj ~ 8.129 x 10^"^, 
a ~ 6.050, ci ~ -2.862, and C2 ~ 2.393. Given S,ciitci), 



we obatin ^(t) directly from Eq.(l) for all the values of t 
where our quantum Monte Carlo data are available. 

In Fig. (1) we present our results. We see a good 
agreement with the theory for ^(i) ^ 6. However, our 
data clearly deviate from the theory for larger values of 
^. In fact the theory overestimates the actual numerical 
data more than 500 percent at t ~ 0.21. We thus con- 
clude that the validity of the new approach is limited to 
very high temperature regime, i.e., T/J > 0.4 for the 2D 
S = 1/2 QHA. 

We would like to thank A. Cuccoli, V. Tognetti, R. 
Vaia and P. Verrucchi for helpful comments and for send- 
ing us their numerical estimates for the function 9'^. 
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Figure Caption: Comparison of the Monte Carlo 
data with the theoretical prediction Eq.(l) for the 2D 
S — 1/2 QHA. The statistical errors of our data are 
much smaller than the size of the symbol. 



[1] A. Cuccoli, V. Tognetti, R. Vaia, and P. Verrucchi, Phys. 

Rev. Lett. 77, 3439 (1996) 
[2] S. Chakravarty et al, Phys. Rev. B 39, 2344 (1989) 
[3] B. B. Beard and U.-J. Wiese, Phys. Rev. Lett. 77, 5130 

(1996) 
[4] J. Apostolakis et al, Phys. Rev. D 43, 2687 (1991); J.-K. 

Kim, Phys. Rev. D 50, 4663 (1994); S. Caracciolo et al., 

Phys. Rev. Lett. 75, 1891 (1995) 
[5] M. S. Makivi'c and H.-Q. Ding, Phys. Rev. B 43, 3562 

(1991) 



100 



iJLP 



10 



0.20 



theory (Ref.(1)) 
our MC data 
□ MC data from Ref.(4 



0.30 



0.40 



